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ABSTRACT 

We report measurements of the X-ray expansion of the youngest Galactic supernova remnant, Gl.9+0.3, using 
Chandra observations in 2007, 2009, and 2011. The measured rates strongly deviate from uniform expansion, 
decreasing radially by about 60% along the X-ray bright SE-NW axis from 0.84% zb 0.06% yr -1 to 0.52% d= 

0.03% yr -1 . This corresponds to undecelerated ages of 120-190 yr, confirming the young age of Gl.9+0.3 and 
implying a significant deceleration of the blast wave. The synchrotron-dominated X-ray emission brightens at a rate 
of 1.9% =b 0.4% yr -1 . We identify bright outer and inner rims with the blast wave and reverse shock, respectively. 
Sharp density gradients in either the ejecta or ambient medium are required to produce the sudden deceleration of 
the reverse shock or the blast wave implied by the large spread in expansion ages. The blast wave could have been 
decelerated recently by an encounter with a modest density discontinuity in the ambient medium, such as may be 
found at a wind termination shock, requiring strong mass loss in the progenitor. Alternatively, the reverse shock 
might have encountered an order-of-magnitude density discontinuity within the ejecta, such as may be found in 
pulsating delay ed-detonation Type la models. We demonstrate that the blast wave is much more decelerated than 
the reverse shock in these models for remnants at ages similar to Gl.9+0.3. Similar effects may also be produced 
by dense shells possibly associated with high-velocity features in Type la spectra. Accounting for the asymmetry 
of Gl.9+0.3 will require more realistic three-dimensional Type la models. 
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1. INTRODUCTION 

Gl.9+0.3 is the remnant of the most recent supernova (SN) 
in our Galaxy (Reynolds et al. 2008, Paper I). Its expansion of 
about 1.6% between 2007 and 2009 measured with Chandra 
(Carlton et al. 2011, Paper V) gives an expansion (undeceler- 
ated) age of about 160 yr, but the estimated mean expansion 
index m (R a t m ) of about m = 0.7 gives an age of ~ 100 yr. 
The X-ray spectrum is dominantly synchrotron with high ab- 
sorption (Reynolds et al. 2009, Paper III); lines of Fe, Si, and 
S are found in small regions with spectroscopic velocities of 
about 14,000 km s -1 (Borkowski et al. 2010, Paper IV), consis- 
tent with the proper motions for a distance of the order of 10 kpc 
(Roy & Pal 2014). For an assumed location near the Galactic 
Center (< d = 8.5 kpc), the diameter is about 2 pc. Gl.9+0.3 is 
quite asymmetric (Figure 1), and the thermal emission is also 
asymmetrically distributed (Borkowski et al. 2013b, Paper VI). 
The relative prominence of the Fe lines, high velocities, the 
absence of a pulsar- wind nebula, and the bilaterally symmetric 
synchrotron emission (as in SN 1006) all point toward a Type la 
origin. Only a very unusual core-collapse event could reproduce 
the observations, while a reasonable SN la model can reach the 
observed size and velocity for a mean external density of about 
0.02 cm -3 (Paper V). Observing the rapid evolution of Gl .9+0.3 
in morphology and brightness can provide unprecedented infor- 
mation concerning the dynamics of SN ejecta and on particle 
acceleration. 

2. OBSERVATIONS 

Chandra observed Gl.9+0.3 with the ACIS S3 chip over 

three epochs: (I) 2007 February and March, (II) 2009 July, and 


(III) 2011 May and July (details can be found in Papers I and 
III- VI). All observations were done in Very Faint Mode, and 
reprocessed with CIAO version 4.6 and CALDB 4.5.9. The 
corresponding effective exposure times are 49.6 ks, 237 ks, 
and 977 ks. Exposure- weighted time intervals between the deep 
Epoch III and shorter Epoch I and II observations are 4.274 yr 
and 1.861 yr. Alignment of observations at different epochs is 
performed simultaneously with expansion measurements. 

Images 512 2 pixels in size were extracted from the merged 
event files by binning event positions to half the ACIS pixel 
size, so one image pixel is Of ! 246 x (X/246. We also extracted 
data cubes 512 2 x 64 in size using the same spatial pixel size, 
while spectral channels from 84 to 595 (1.2-8. 7 keV energy 
range) were binned by a factor of eight. 

X-ray spectra were extracted from individual rather than 
merged event files and then summed (response files were av- 
eraged). Spectral analysis was performed with XSPEC version 
12.8.1 (Arnaud 1996), using C-statistics (Cash 1979). The back- 
ground was modeled rather than subtracted. Spectra of Gl.9+0.3 
were modeled with an absorbed power law, using the solar abun- 
dances of Grevesse & Sauval (1998) in the phabs absorption 
model. Thermal emission contributes negligibly to broadband 
fluxes. 

3. EXPANSION 

We first measured the overall expansion using the method 
described in Paper V. Briefly, we smoothed the 2011 Chandra 
data cube with the spectro- spatial method of Krishnamurthy 
et al. (2010), and summed the smoothed cubes in the spectral 
dimension to arrive at smoothed images in the 1.2-8 keV 
energy range. (Softer X-rays are absorbed by the intervening 
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Table 1 

Expansion Rates and Fluxes 


Region 

Expansion Rate 
(% yr -1 ) 

S 20 07 a 

(% yr 

S2009 a 

N H b 

(10 22 cm -2 ) 

r c 

F 2007 d 

(10 

E 2 009 d 

-13 —2 - 

J erg cm ^ s 

f2011 d 

‘) 

Flux Rate Increase 
(% yr -1 ) 

Total 

0.589 

-0.2 

0.4 

7.25 

2.40 

27.34 

28.79 

29.73 

1.9 


(0.573, 0.605) 

(-0.7, 0.2) 

(-0.1, 0.9) 

(7.16,7.34) 

(2.37, 2.43) 

(26.75, 27.93) 

(28.50, 29.08) 

(29.57, 29.99) 

(1.5, 2.3) 

Outer SE-NW pair 

0.523 

0.0 

1.9 

7.90 

2.13 

2.72 

2.87 

2.89 

0.9 


(0.487, 0.560) 

(-1.5, 1.7) 

(0.2, 3.6) 

(7.60, 8.21) 

(2.04, 2.23) 

(2.61,2.83) 

(2.81,2.92) 

(2.86, 2.91) 

(-0.3, 2.2) 

Middle SE-NW pair 

0.616 

-0.4 

1.0 

7.58 

2.34 

9.16 

9.75 

10.05 

1.9 


(0.592, 0.640) 

(-1.1, 0.4) 

(0.2, 1.8) 

(7.45, 7.72) 

(2.29, 2.38) 

(8.97, 9.35) 

(9.66, 9.84) 

(10.00, 10.09) 

(1.3, 2.5) 

Inner SE-NW pair 

0.842 

0.3 

-0.3 

7.02 

2.38 

3.56 

3.66 

3.91 

2.8 


(0.783, 0.898) 

(-0.9, 1.6) 

(-1.6, 1.0) 

(6.82, 7.23) 

(2.31,2.45) 

(3.45, 3.70) 

(3.60, 3.71) 

(3.88, 3.94) 

(1.8, 3.8) 

N-S pair 

0.576 

-0.2 

-0.2 

7.09 

2.51 

8.23 

8.65 

9.00 

2.1 


(0.544, 0.609) 

(-1.0, 0.7) 

(-1.1, 0.6) 

(6.96, 7.23) 

(2.46, 2.56) 

(8.05, 8.41) 

(8.56, 8.74) 

(8.95, 9.05) 

(1.4, 2.8) 


Notes. Expansion rates and fluxes in odd rows, 90% confidence limits in even rows. 
a Surface brightness change. 
b Hydrogen column density. 
c Power-law photon index. 
d Absorbed flux in the 1-7 keV energy range. 


interstellar medium, while the background dominates at the 
highest energies.) We used the smoothed image as a model 
for the surface brightness of Gl.9+0.3 at Epoch III (Figure 1). 
This model image was background- subtracted and then fit to the 
unsmoothed 1.2-8 keV images from earlier epochs (i.e., shrunk 
to fit) using C-statistics. Seven point sources within Gl.9+0.3 
were masked out. There are four free parameters in this model: a 
physical scaling factor, a surface-brightness scaling factor, and 
expansion center coordinates. Independent fits to Epoch I and 
II images were consistent with the constant expansion. We then 
assumed the same expansion rate while fitting jointly to the 
Epoch I and II images, but allowed for independent surface- 
brightness scaling factors and expansion center coordinates. 
The results of all of our measurements are listed in Table 1. 
The measured expansion rate is 0.589% ±0.016% yr -1 (all 
errors are at 90% confidence), in reasonable agreement with our 
previous measurement of 0.642% ± 0.049% (Paper V). Within 
the errors, the mean surface brightness does not vary with time, 
again in agreement with Paper V. 

To search for deviations from uniform expansion, we system- 
atically examined the motions of large, distinct spatial struc- 
tures. The final regions chosen are overlaid atop X-ray and 
radio images in Figure 1 . They include a pair on the N and S 
rim and three pairs of regions along the major SE-NW axis: an 
outer pair (“ears”), a middle pair containing the brightest rim 
emission, and an inner pair with distinct interior rims. 

Measurements of the expansion for each region pair were 
performed in the same way as for the entire remnant. We 
assumed a common expansion center of coordinates for all 
of the pairs for Epochs I and II, but the measured expansion 
rates are insensitive to this assumption. Slightly different region 
positions and areas (but not shapes) were used for Epochs I 
and II, to approximately match the measured expansion rates. 
The common expansion centers of the (J2000) coordinates are 
R.A. 17 h 48 m 45^639 ± 0 S .004 (17 h 48 m 45 s .629 ± 0 S .002), decl. 
— 27°10 / 06 / . / 85 ± 0?06 (-27°10'06?96 ± 0T03) for Epoch I (II). 
A small (Of! 11) but significant difference between these centers 
suggests that the coordinates of the corresponding reference 
observations ID6708 and ID10112 are slightly misaligned. The 
magnitude of this shift is consistent with the Chandra external 
astrometric errors (mean error of Of! 16; Rots 2009). 

The measured expansion rates strongly deviate from uni- 
form expansion (Table 1). Expansion rates increase inward 


by about 60% along the bright SE-NW axis, ranging from 

0. 52% ± 0.03% yr -1 for the outer ears to 0.84% ± 0.06% yr -1 
for the inner rims. The bright rims in the middle expand slightly 
faster (0.616% ± 0.024% yr -1 ) than the ears, but even this small 
difference is statistically highly significant. The brightness- 
weighted linear displacement is Of! 29 yr -1 for all three rims. The 
N-S expansion is intermediate (0.58% ± 0.04% yr -1 ) between 
the bright middle rims and the ears, with an average displace- 
ment of Of! 23 yr -1 . We demonstrate this differential expansion 
for two representative profiles in the SE and NW based on the 
Epoch II observations (Figures 2 and 3). As can be seen most 
clearly in the close-up inserts, the expansion rate (green) that 
matches the bright middle rims is too slow (too little shrinkage) 
for the inner rims, while the faster expansion rate (red) which 
matches the inner rims is too fast for the middle and outer rims. 

Contributions to the C- statistic along the profiles quantify 
the fit quality. Systematic deviations are present even for the 
best-fit models, perhaps due to spatial variation of the Chandra 
point-spread function, smaller-scale spatial variations in the 
expansion rate, or projection effects. Strong counting noise 
limits the accuracy of the expansion measurements. However, 
these profiles are merely illustrative; the results in Table 1 
are derived from the expansion of the entire regions shown 
in Figure 1 . 

4. FLUX INCREASE 

The spatially integrated spectra of Gl.9+0.3 from each epoch 
were jointly fit, together with their background spectra, with an 
absorbed power law in the 1-9 keV energy range. The power- 
law index and absorbing column density were assumed to be 
constant in time, but the fluxes at each epoch were left free. 
There is good agreement with previous measurements for Epoch 

1, but the newly determined Epoch II flux is 2.5% larger, which 
disagrees with the previous measurement also reported in Paper 
V, presumably due to updates in the ACIS S3 calibration that 
were applied to all three data sets. 

Spatially integrated fluxes are clearly increasing with time 
(Table 1). A likelihood ratio test reveals that the linear flux 
increase is consistent with the individual flux measurements. 
The measured flux increase is 1.9% ± 0.4% yr -1 , in agreement 
with the previous, less accurate value of 1.7% ± 1.0% yr -1 
(Paper V). 
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Figure 1. Top: 2011 Chandra image of G1 .9+0.3. The scale is in counts per ACIS 
pixel in the 1.2-8 keV energy range (image was smoothed with the multiscale 
partitioning method of Krishnamurthy et al. 2010). Bottom: 1.4 GHz Very 
Large Array radio image from 2008 December. The scale is in Jy beam -1 . The 
resolution is 273 x 174. N is up and E is to the left. Intensities are shown with the 
cubehelix color scheme of Green (201 1). Pairs of regions chosen for expansion 
studies are overlaid: outer (green), middle (magenta), and inner (yellow) pairs 
along the SE-NW axis, and NS (red) pair. The image size is 123" x 123". 

We also measured the fluxes and rates of flux increase 
for the region pairs shown in Figure 1 in the same way as 
for the spatially integrated fluxes, except that the background 
contribution was scaled down (rather than fit again) by the 
region/total area ratios from the global fit. For each region pair, 
the measured rate of flux increase is consistent with the spatially 
integrated rate. There is an apparent trend in the brightening rate 



Figure 2. Top: SE profile from 2009 (blue) along the position shown on the 
2009 image (left inset), together with model profiles corresponding to best-fit 
expansion rates of 0.62% yr -1 (green) and 0.84% yr -1 (red) for the middle and 
inner rims (see Figure 1 and Table 1). Right inset: zoomed view of the inner 
rim. Bottom: contribution to C-statistic. The horizontal scale is distance along 
the profile, measured from E to W. 

with radius, but uncertainties are large due to uncertainties in 
expanding the regions and to the dust scattering of emission from 
brighter regions (Paper III), as well as to large measurement 
errors. We conclude that evidence for spatial variations in the 
rate of flux increase is weak, although it should be the subject 
of future investigations. 

5. DISCUSSION 

The nonuniform expansion we observe for our three pairs of 
regions can be rephrased as large differences in the expansion 
age (f exp = At R/AR ), in the sense that f exp is largest (greatest 
deceleration) for the slower-expanding outermost material (see 
Figure 1 and Table 1). The outermost ears, the bright rim 
emission in the middle, and the distinct interior rims all have 
measurably different expansion ages: 190 yr, 160 yr, and 120 yr. 
We define the expansion indices m = d \nr/d \nt so that a 
feature at radius r obeys r a t m (note that d Inr/d \nt = 
d In R/d In t, where R is the projection of the true radius r onto 
the plane of the sky). Then, the true remnant age t = mf w t Qxp , 
where mf w is the forward- shock expansion index. Since we only 
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Figure 3. NW profile from 2009. See the caption of Figure 2 for explanation. 


have an estimate of 100 yr for the true remnant age (Paper V), 
we can only determine relative m values. 

The large spread in expansion ages between the inner rims 
and the ears implies a large deceleration of the forward shock, 
m fw < 0.6, significantly stronger than expected (mf w ~ 0.7) 
in the models with smoothly varying density distributions for 
ejecta and ambient medium considered in Paper V. However, 
the more fundamental problem is that, contrary to observations, 
m in these models typically varies only slightly (if at all) from 
the reverse shock to the blast wave, and it increases instead of 
decreasing with radius (Chevalier 1982; Dwarkadas & Chevalier 
1998, hereafter DC98). Either the ejecta or the ambient medium 
density distribution (or both) must be very different from the 
slowly varying density distributions considered so far. Various 
possibilities include substantial clumping or sudden density 
jumps within the ejecta or the ambient medium. 

Fine- scale clumping within the ejecta is unlikely to explain 
the faster than expected expansion of the inner rims, however. 
Much like the outer rims, the inner rims consist of cohesive 
and continuous filaments in contrast to the much more clumpy 
north rim (Figure 1) where the ejecta emit most strongly in 
thermal X-rays (Paper VI). Another explanation involves a rapid 
deceleration of the blast wave following a recent encounter with 
a moderate density jump in the ambient medium. In order to 
account for the approximately 60%-70% slower expansion of 
the middle and outer rims relative to the inner rims (Table 1), 


a factor of 3-6 density jump is required (from Equations (4) 
and (6) in Borkowski et al. 1997). A density jump of this 
magnitude is consistent with a wind termination shock. If this 
were true, then the SN progenitor must have been losing mass 
in a strong and asymmetric stellar wind. We estimate the wind 
parameter (= M/ v w , where M is the mass-loss rate and v w 
is the wind speed) at 4 x 10 -6 M o yr -1 /l(X)0 km s -1 , using 
the analytic wind thin- shell solution for an exponential ejecta 
density distribution from Paper V with the remnant’s radius of 
2 pc and undecelerated age of 120 yr, and assuming a standard 
thermonuclear explosion with a kinetic energy of 10 51 erg and 
an ejected mass equal to the Chandrasekhar mass. Such a strong 
wind favors a single-degenerate progenitor (Hachisu et al. 1 996). 
The deceleration is m = 0.88 in this model, suggesting that the 
explosion occurred sometime in the first decade of the 20th 
century. In this scenario, it is difficult to understand the origin 
of the strong north-south asymmetry seen at radio wavelengths 
(Figure 1) and in the spatial distribution of the thermal X-ray 
emission. As discussed in Paper VI, a strongly asymmetric Type 
la explosion provides the best explanation for this asymmetry. 

Unlike the case for smooth ejecta density profiles (such 
as those resulting from delayed-detonation explosions), it is 
possible to obtain less deceleration inward (mf w < m rev ) in one- 
dimensional (ID) numerical models for SNe la expanding into 
a uniform ambient medium for ejecta profiles with substantial 
structure (e.g., the PDDe model in Badenes et al. 2003). 
Prominent density structures in the outer ejecta layers are present 
in the deflagration models (e.g., the W7 model of Nomoto et al. 
1984), sub-Chandrasekhar explosions (DC98; Badenes et al. 
2003), and in pulsating delayed-detonation (PDD) models (e.g., 
Dessart et al. 2014, see their PDDEF1 model plotted in Figure 4). 
This density profile can be satisfactorily approximated by a 
power law p oc r~ n with n — 5.8 above the transition velocity 
v tr = 14,800 km s -1 , and by the exponential model of DC98 at 
lower velocities. At u tr , there is a large (1 1.4) density jump. 

In order to investigate how a density jump affects the 
speeds of the reverse and forward shocks, we performed ID 
hydrodynamical simulations using the VH-1 hydrocode (see 
Warren & Blondin 2013 for a recent description of this code) 
for this composite power-law-exponential ejecta model. Here, 
we discuss the results for a preshock density no of 0.2 cm -3 , 
matching the preshock density found in a young Type la SNR 
0509-67.5 in the Farge Magellanic Cloud (Williams et al. 2011), 
but they can be scaled to any value of no (DC98). Prior to the 
impact of the reverse shock with the density jump, the density 
profile can be described by the self- similar solutions of Chevalier 
(1982) with mf w = m rev = (n — 3)/ n = 0.48 (see the density 
profile at 50 yr in Figure 4). The reverse shock arrives at the 
density jump at t = 83(fto/0.2 cm -3 ) -1 / 3 yr with a velocity 
of v s = 3v tr /n = 0.52% (in the frame of reference moving 
with the ejecta), and then splits into transmitted and reflected 
shocks. The transmitted shock is the decelerated reverse shock 
(inward-facing). Its velocity is v t = (f$/8) l/2 v s , where /3 is the 
pressure enhancement that depends only on the density jump 8 
and varies between 1 and 6 (see Equation (6) in Borkowski et al. 
1997). In the rest frame of the explosion, the deceleration of 
the transmitted reverse shock becomes m rev = 1 — 3 (/3/8) l/2 /n. 
With 8 = 11.4, the overpressure is /3 = 2.76 and m rev increases 
to 0.75, in good agreement with hydrodynamical simulations 
(Figure 4). The reflected shock propagates first back into the 
shocked ejecta and then into the shocked ambient medium. At 
100 yr, it has already passed through the low-density contact 
discontinuity that separates the shocked ejecta from the shocked 
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Figure 4. Top: density vs. radius in ID hydrodynamical simulations with 
a composite power-law-exponential ejecta model and the uniform ambient 
medium with no = 0.2 cm -3 at 50 yr (solid blue line), 100 yr (green), and 
150 yr (magenta). Undecelerated (freely expanding) ejecta at 50 yr are also 
shown: the composite model (dashed line) and the PDDEL1 model (Dessart 
et al. 2014) ( + signs). Bottom: radii of the reverse and forward shocks (in blue 
and green solid lines), and their deceleration parameters m (dashed lines). The 
reverse shock expansion is faster than the blast wave expansion during a long 
(~ 100 yr) period of time. 

ambient medium, and can be seen in Figure 4 as a small density 
discontinuity near the contact discontinuity. The reflected shock 
strengthens with time (see density profile at 150 yr), and 
eventually merges with the blast wave at 175 yr, resulting in 
an abrupt increase of mf w from 0.48 to 0.62. By this time, m rev 
has already decreased to about the same value after a transient 
phase following the sudden deceleration of the reverse shock. 
This decrease continues in the subsequent evolution, while nif w 
remains about constant, so that mf w > m rev . 

During an extended period of time comparable to a young 
Type la remnant’s age, PDD explosions show a more rapid 
expansion of the reverse shock than the forward shock 
(Figure 4). While this still might be a reasonable interpretation 
for Gl. 9+0.3, published PDD models fail to match its proper- 
ties in detail. Assuming the ejecta structure of the composite 
model just described and identifying the blast wave with the 
bright middle rims at r = 2.04(d/8.5 kpc) pc expanding at 
0.62% yr -1 , the inferred preshock density no and free expan- 
sion ejecta velocity at the reverse shock are 0.08(d/8.5 kpc) -5,8 
cm -3 and 18,000(d/8.5 kpc) km s -1 . At 8.5 kpc, no is several 
times higher than estimated in Paper V, but the ejecta velocity 
is consistent with previous estimates. Gl.9+0.3 is dynamically 
too young, however, for the reverse shock to have reached the 
density jump at u tr = 14,800 km s -1 . No published PDD model 
has Vtr as high as 18,000 km s -1 , with the highest (16,000 km 
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s -1 ) found in the PDDa model (Badenes et al. 2003). Among 
the models considered by Badenes et al. (2003), only the sub- 
Chandrasekhar model SCH has a higher v [r at 18,500 km s _1 . 
Furthermore, sub-Chandrasekhar explosions produce large (up 
to several x 10 _3 M o ) amounts of radioactive 44 Ti (e.g., Woosley 
& Kasen 2011), while the strength of the (possible) 44 Sc line 
implies that at most ~ 10 _5 M Q of 44 Ti was expelled by the SN 
that produced Gl.9+0.3 (Borkowski et al. 2013a). 

Only the outermost ( v > 18,000(d/8.5 kpc) km s -1 ) ejecta 
have been shocked so far in Gl.9+0.3 (velocities might be 
as high as v > 21,000 km s -1 if d > 10 kpc; Roy & Pal 
2014). Discrete absorption features with such extreme velocities 
are commonly seen in Type la SN spectra (e.g., Childress 
et al. 2014), perhaps indicating the presence of dense shells 
of material (Tanaka et al. 2006). Their origin is unknown, 
but possible explanations include a collision of ejecta with 
circumstellar medium or clumping within the ejecta. Fast dense 
shells with an order- of-magnitude density contrast relative to 
overlying tenuous ejecta are expected to decelerate the reverse 
shock just like the density jump in PDD models. This might be 
the explanation for the fast expansion of the reverse shock in 
Gl.9+0.3. 

The flux increase may be due to global effects such as an 
increase in magnetic-field strength, which would raise radio 
fluxes as well, and/or an increase in the maximum energy to 
which electrons are accelerated, which could raise X-ray fluxes 
even with constant radio. The radio flux of Gl.9+0.3 is rising at 
l%-2% yr -1 (Green et al. 2008, Paper II; Murphy et al. 2008); 
our ongoing JVLA observations should reduce the uncertainty 
involved in determining whether this rate is consistent with 
the X-ray rate. A higher X-ray than radio rate of increase 
would require an increasing maximum energy. This would be 
a very interesting result. Spatial variations contain additional 
information. One might naively expect the greatest brightness 
increase at the location of greatest deceleration (thermalization 
of kinetic energy); our data suggest but do not compel this 
conclusion. 

Gl .9+0.3 offers us a unique opportunity to study the dynamics 
and spatial distribution of the outermost ejecta of a likely Type 
la SN, and the process of particle acceleration in very fast 
collisionless shocks. Our current understanding of its dynamics 
and of the spatial variations in the flux rate increase is limited 
by the short-time baseline of the Epoch I— III observations and 
the poor signal to noise of the 2007 and 2009 Chandra data 
sets. Further Chandra observations will allow for significant 
advances in our understanding of Gl.9+0.3. 

This work was supported by NASA through Chandra General 
Observer Program grants SAO G01-12098A and B. 
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